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A NEW METHOD OF COMPUTING PENETRATION VARIABLES FOR 
SHAPED=CHARGE JETS 


In earlier attempts to describe the penetration of shaped-charge jets, 


the total depth of penetration was computed from J vat. Evaluation of the 
integral has not been possible because the penetration velocity U for the 
discontinuous portion of the jet has not been properly defined. In the 
following analysis, it is postulated that the total penetration can be 
represented by { vat +S ap, where the integral is taken over the penetration 
produced by the continuous portion of the jet and the summation is taken over 
the penetration produced by the discrete particles resulting from jet breakup. 


It is shown that numerical values of the penetration variables for the 105mm 
unconfined test charge are compatible with the modified penetration equation. 


Scaling relations are introduced, and the penetration equations are given in 


a 
it was assumed that the incompressible approximation provides an adequate 


relationship between the jet and penetration velocities for a continuous jet, 


4+ dai ila ahnum that tha Ff 
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materials do not seriously alter this relation. The incompressible approximation 


does underestimate the pressure by an appreciable amount. 
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INTRODUCTION 


A shaped charge is a jet-forming explosive device containing an axially 
symmetric cavity lined with a thin metal liner. The geometrical design of a 
shaped charge can be varied over a wide range; however, cylindrical charges 
containing conical cavities have been most extensively investigated. Line 
drawings illustrating the jet formation end penetration process for a typical 
shaped charge are presented in Figure 1. When the explosive is detonated, the 
high pressure of the detonation products collapses the metal liner upon itself 
forming a high-velocity, forward-moving, metallic jet, which is able to 
penetrate high-strength materials such as erin: Different elements of the 
liner experience different collapse velocities and the resulting jet contains 


with 


with forward portions of the jet traveling faster 
eae) 


a spectrum of velocities with forward portions 


than succeeding portions. The velocity gradient in the jet causes it to 
lengthen as it travels forward from the charge. Initially the jet lengthens 
in a ductile manner, but eventually the ductility of the material is exceeded 


and the jet fractures into a series of discrete particles. 


The pressure exerted by the impinging jet greatly exceeds the yield 
strength of the target material, and penetration takes place by a hydrodynamic 
process as illustrated in Figure 1. Theories describing the penetration of a 
shaped-charge jet were advanced independently by Pugh? in the United States and 
by Hill, Mott, and Pack” in Great Britain. Both theories were based on steady- 


state incompressible fluid flow and resulted in an equation of the form 


ny 


no ,(V, - ij>-= o,U", (1) 


where 0. is the jet density; Py» the target density; Vag the velocity of the 
jet at the stagnation point; U = ies the penetration velocity; and i, a 
in the transfer of momentum 


to the target from a continuous jet and from a jet of discrete particles. The 


accuracy of the penetration equation has been extensively investigated by 


2K 


The existence of independent work by Hill, Mott, and Pack was cited in 


picheibeszer’, who further generalized (1) by including an additive strength 
term. Attempts to use (1), or some modification of it, to correlate experimental 


values of U and ve have been extremely successful. 


By introducing plausible assumptions concerning the behavior of the jet in 
flight, Pugh was able to integrate the penetration velocity and account for 
qualitative features of the penetration-standoff curve. This method requires 
defining an appropriate penetration velocity U for the penetration produced by 
the jet after it has broken into a series of discrete particles. Because the 
correct definition of U is not immediately obvious, this report describes a 
different approach to the problem in which the penetration process is 
represented by different mathematical models depending on whether the jet is 


continuous or broken-up. 
AN APPROXIMATE SOLUTION FOR THE PENETRATION VARIABLES 


To calculate the total depth of penetration in @ massive target, it is 
necessary to add the penetrations produced by the various jet elements. For 
a@ continuous jet, correct values are obtained by integrating the penetration 
equation. After the jet fractures into a series of discrete particles, the 
physical process is better regarded as an addition of discrete elements of 
penetration corresponding to discrete particles of jet, and a more appropriate 
mathematical description is a summation of finite increments. Accordingly, 


the total penetration P_ is given by 


T 
ty N 
Pe -[ Udt + ) AP, , (2) 
t isl 
1?) 


where ve is the time at which the jet reaches the target; ti» the time at which 
the jet breaks up; and AP. 5 the penetration produced by the ih particle. 


For a standoff of two or more cone diameters, it is well know that the 
velocity of the forward tip of the jet emerging from a target of thickness 
P,(0< P< Pa)» can be calculated (to within a few per cent) from penetration 
time measurements by assuming that all jet elements are formed simultaneously 
at a virtual origin. Therefore an approximate solution for the penetration 


variables can be obtained by assuming that, at t = 0, the jet is formed at a 


8 


virtual origin located a distance z_ in front of the target surface. As 


°o 


nf «4 ot os Pega oe <r aon co a sro laAnmtekar ama 
shown in Figure 2(a), the penetration depth and impinging jet velocity are 


related by 
(3) 


where v(t) is the velocity of the jet at the stagnation point. It is further 


assumed that the effects of strength are negligible; hence, for a continuous 


wet gOS 


jet (x = 1), Equation (1) reduces to 
U=Vv,/(1 + 7), (4) 

i yi/2 ow fw e ~~ lat a Los ay. ae Se ers me 
where y = (p p,/0,) "6 hen (3) is differentiated with respect to time, one 
obtains 

U = dP/dt = Vv, + tav ,/dt. (5) 
After the penetration velocity is eliminated between (4) and (5), the resulting 
equation in V; and t can be integrated to give 
° i y/(1 7 + \ perk 
se (oe ie (6) 


ich it first encounters the target, relative to t = 0 at the virtual 


1 %f¢7N 5 Ee eae 


: . — ss am Boke gabe a ae. dete ete mm nesse m e wede  A 
origin. Equation (6), together with (3), can be used to compute penetration- 
o 
e 


determined experimentally, where Equation (3) is now given oy. 
Oo 1 
P(t) = V9 & (t/o) 7/0 +9) (1 


Af 5 
ALE a 


continuous jet of discrete particles. An experimental 


O ts 


jet 
determination of t, can be obtained from the point at which the penetration- 
aS 


urve departs from that predicted for a con 


tinuous jet. If it is assumed 


that break-up occurs simultaneously throughout the remaining jet, then, as 
shown in Figure 2(b), the following relation must hold: 


\o 
a 


DE. es EE Ta (al! 


mim 
doe 
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material and a is the length of the effective jet that has not yet 
=1: 

penetrated the target. To calculate the total penetration, it is assumed that 
f 


ea 


yp 
re) 
x 
p 
4 
j--> 
Q 
+ 
(D 
O 


i.e., AP. = AL. |. According to (2) and (8) and the definition of AP, 5 the 
total penetration will be given by 


~— been | 
= : in 
P, = P(t,) + E + P(t,) - +07" | /7; (9) 
which is independent of the total number of jet particles involved. 
Flash radiographs indicate that the broken jet can be regarded as a series 
of discrete particles of approximately equal ieweth?. If it is assumed that 
AL, = AL,= cee Aly = Af, the number of effective jet 


N = [z., + P(t.) - evn | /Qs. (10) 


Fiash radiographs also show that the increment in velocity between adjacent 
= 
particles is approximately a constant~. Hence, the increment in velocity 


between adjacent particles is given by 


W = - ly (t+_) - ymin | /(w -~ 1) (11) 
J | J° 1’ J | 7oN\* banat A Van] 
A th 7 + f +h th +4 h 
ana Uae VeiocL Uy OL Loe L parvicie LS given DY 
Vis = Vialttsy: + Ge dav... (12) 
Jol Jot d 
An inspection of Figure 2{c) shows that the time between the arrival of 
i .th = in a 74 s a a ,st ij 2 bo} e s 7. 
the i” and the (i +1) particle is given by 
Abs’ a = |AP + Ag - AV.(t, - t ie ee He 
i,i+l $ i 1) /( pin (15) 


Of course, the arrival time of the oe particle is given by 


i 
t, = +, + ) At ay > (14) 
— : 


4-2 
JT 


where i>e. 


The penetration depth at the instant the a particle reaches the bottom 
of the hole is given by 


P(t, ) = P(t, ) + (i - L)AP. (15) 


Tear 


The penetration by the discontinuous portion of the jet will be completely 
specified once the constants Aé# and i have been determined. The number of 
particles is determined from (10), and their respective velocities are 
determined from (11) and (12). The time at which each particle starts to 
penetrate the target can be determined from (13) and (14), and the corresponding 
penetration depths can be obtained from (15). 


The accuracy of the theory for correlating the penetration variables is 
best illustrated by analyzing data obtained with the 105 mm unconfined test 
charge used at the Ballistic Research Laboratories. The time at which fracture 
of the jet occurs was determined by integrating the hydrodynamic penetration 
equation for a continuous jet and comparing it with the experimental penetration- 
time date . The point at which the two diverge; i.e., the time when the 
continuous jet approximation no longer adequately describes the penetration 
process, is taken to be the time at which simultaneous fracture occurs in the 
remainder of the effective jet. As shown in Figure 3, the divergence of the 
experimental data from the continuous jet approximation occurs at 120 u-sec. 
The value of AZ was estimated from flash radiographs to be 1.68 em. The values 
of ce and v; were measured directly and found to be 0.210 em/p-see and 


0.701 em/y-sec, respectively. 


The penetration theory for discrete particles can be used to extend the 
penetration-time curve beyond 120 p-sec. As shown in Figure 3, results of the 
calculation are in good agreement with experimental observations. In addition 
to the penetration-time data, the velocity of jets emerging from targets of 
various thicknesses gives an independent measure of the jet velocity as a 
function of the penetration depth. As illustrated in Figure 4, the curve 
computed from the penetration theory adequately represents the emerging jet 


velocity measurements. 


See reference 4, for a general description of the experimental techniques 
used in obtaining penetration-time data. 


LL 


THE USE OF CONTINUOUS VARIABLES TO REPRESENT 
THE PENETRATION PRODUCED BY A PARTICLE JET 
It has already been noted that the total penetration produced by the 
broken-up portion of the jet depends on the total length of all the particles 
involved, but not on the length of an individual particle. It is therefore 
permissible to represent the penetration variables by continuous functions 
obtained from appropriate mathematical limits taken as the particle length AZ 
approaches zero and the number of particles N approaches infinity. Of course, 
the limits must be taken in such a way that 


L= Lim (NAL) = 2. + P(t,) - t1V (16) 
As» 0 ° 1 a oa 
N—@ 


is a constant. 


An expression for the velocity of the impinging jet is obtained after 


dividing (11) by Ag, 


Re Dnt ig as pe min) ,, 
Napoe ie [' 38°) re ve | / (NAL re Ag); (17) 
taking the limit as A# tends to zero, 
Lim (Av ,/AZ) = av./dé = - ly yin | /L (18) 
Als O J J Ld 4 J -3 
and integrating, 
Vs) s 
lan rT ; Th gy 
av =~ (1/L) vat) - "| dé. (19) + 
3 fy) 0 


The final expression for the jet velocity is 


v.(s) = Vette) - Pa ye ynin| s/L (30) 
J. é qo eed [ ae 1? j a / 3 { ) 
where the rameter s = L £is the location of an infinitesimal 
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An expression for the penetration is readily obtained by taking the limit 


of AP, as AZ tends to zero, 
Lim (Ap/Ag) = aP/ae = 1/y (21) 
Ag +0 
and integrating, 
P(s) la 
| aP = Q/y) | af. (22) 
P(t, ) @) 
The resulting expression for the penetration is 
P(s) = P(t,) + 8/7. (23) 
In much the same way, an expression for the time is obtained after 
dividing (13) by Ag. 
(V5) yp (Oty 34/44) = (AP/AL&) +1 - (A V ,/O8) (t, - t)) » (24) 
taking the limit as A# tends to zero, 
V ,at/as = (1/7) +1 - (t, - t, av ,/as ; (25) 
and integrating 
v ,(s)t(s) 8 V ,(s) 
{ a(v,t) = (1+y7) dg +t. av... (26) 
a 7 J a d 
Vee 0 v.(t,) 
a aes mae J 
The final expression for the time is 
t(s) = te s(y + 1)/W ,(s). (27) 


Equations (20), (23), and (27) are parametric equations from which the 
penetration variables can be calculated. The parameter is permitted to vary 
continuously over the range 0 <s <L, and the penetration variables calculated 


from the equations are continuous functions. A graph illustrating the 


functional relation between the penetration and the time is presented in Figure 5. 


For comparison purposes, results computed from the discrete particle theory 


have been superimposed. The time and penetration depth corresponding to each 
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particle impact on the target were computed from (14) and (15) respectively. 
To emphasize the discrete nature of the process, the time required for each 


particle to complete its penetration was computed from 
= = ; 28 
At, = AP, /U, = (1 + 7) P/V), (28) 


According to the idealized model of the particle jet, there will be intervals 

of time (Obs aaa - At, ) during which no penetration takes place The length 

of this interval increases with the depth of penetration because adjacent particles 
have more time to spread due to differences in their velocities. As shown in 
Figure 5, the continuous parameter representation of the penetration-time curve 
follows the discrete particle theory very closely. Because the equations are 
simpler, future calculations of penetration variables will use the continuous 


parameter representation. 


For a continuous jet Bernoulli's equation was used to establish a relation 
between the two velocities V. and ut. A similar relation can be obtained for a 
jet of discrete particles by differentiating (20), (23) and (27) with respect to 
t, obtaining 


av ,/at = + [vjce) - va (ds/at)/L , (29) 
U = ap/at = (1/y) ds/dt , (30) 
(t - t, av ,/at + Vv, = [ + 1)/7] ds /at (31) 


respectively. When (30) is used to eliminate ds/dt from (29) and (31), one 
av ,/at 2g [v (ey) 2 vain] yu/L , (2) 


(es t, av ,/at t= (y + 1)uU. (33) 


The discontinuous nature of the penetration variables was neglected in 
Figure 3 and Figure 4, a smooth curve being drawn through points representing 
the variable at the instant each particle impacts with the target. 
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Elimination of av ./at between (32) and (33) gives 


n [ ,(¢) 7 in| yu(t - t,)/L + v, = (y + 1)uU. (3h) 


To eliminate (t - t,) from the preceding equations, s is eliminated between 
(20) and (27) to give 


| “ V ° 
sh Vamauas Pe ce j | 


Finally, substitution of (35) into (34) gives 
e : Pear 
U= v/(y € 1)v ,(t,). (36) 


For a jet of discrete particles, the penetration velocity is proportional to 


s 2 
the square of the jet velocity, and for t 


h the y, and s reason, the properties of a 


fr 


discontinuous jet cannot be described by replacing O, in Equation (1) w 


effective density. It is interesting to note that the problem of defining an 


iF + 
t Ll 
{t) = J ae a J ars (be 2" Jt) 
G v u 
ro) fe) 1 


Where U, is given by (4) and Up by (36). 


SCALING RELATIONS 


Ky 


t has been shown experimentally that the penetration-time curves for a 
t 


es of scaled charges reduce to a single curve when both the 


penetration and time are divided by the charge diameter’. The accuracy with 
which a single curve represents the penetration-time data is illustrated in 
Figure 6. For the data presented in Figure 6, the location of the virtual 
origin with respect to the top of the target can be expressed as Dis ee ¢ oF? 
where D is the charge diameter and bo? a constant. 


| a | 


bos PART Dee arena fz « fL\% tame =e ey ae 
Equations (3) and (6) can be used to snow 


in charge diameters, is given by 


where 


(x /)/4 + 7) (39) 


and + = t/D. (40) 


After the jet breaks up the penetration can be expressed as a function of the 
dimensionless parameter go = s/D. The penetration, measured in charge diameters, 


is given by 


I 


P/D = x(a) x(t, ) + of/y (41) 


and the corresponding times are given by 


t/D = t(o) = tT, + aly + L/W, ; (42) 
where V. is now given by 
J 
me , = ‘ae min 7 min 
V; =V (J [ vj) V | o/ E + x(t, ) TV, : 
f)zy 
40) 


i an aot, 2 In ana 
lables, 7 = t/D and x = P/D, and 


), (20), (23) and (27). For a homologous 
a 


series of scaled charges there exists 


are obtained from Equations (1 


mis 


unique relation between n and Tt, which 
t 


size. velocity U an 


5 
of 
velocity Va can also be ae as functions of either n or 1, the functions 


being cag ieee of charge size. 
GENERALIZATION OF THE THEORY FOR SHORT STANDOFF 


The present theory will be inaccurate for short standoffs where it is no 

ié to neglect small differences between the times and positions 
at which different jet elements are formed. To remove the inaccuracies at 

short standoffs, a complete description of the linercollapse and jet-formation 
process is required. At the present time the necessary data are available for 
only one charge, the 105 mm unconfined test charge used at the Ballistic 
Research Laboratories’. For this charge, the collapse velocity in em/y-sec 


imparted to an element of the liner is give 


* ¥ he 
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ve =- 0.197898 + 0.288477x- 0.0775152x° +0. 00942615x7 


~ 0.000445 36. *, for 4.0 em <x < 8.0 cm, (44) 


where x is the axial position of the liner element measured relative to the 


apex of the conical liner. 


mee 2 am: 2 = oe ae ee ee on EW PE Tr _—_ a ae ae ad 


€Imsqpntanna Ananwnethine o meant See be ea. At wre lh +o +4 * a 
BQUAtULONS Gesecr1lioings a POLlNG LNLviatea Cnarge OL LeENngsun b CONLAININE a 


e 
conical cavity of height H are readily derived using the quantities illustrated 


“~ Re ieee ye. S0seh eee as 1. ~ 2 cal ee ee ee aS ea ig Sige oe Meas sae eR 
HlO CAPLOUSIVeE Lo Lill ULaALeEU. Lie Li required 1or bu OlL LOW wave O Lear 
a liner element located at the position x is given by 
2 2| 1/2 
T(x) = fe -H+x) + (x tan a) / fees (45) 
_ awk Vv 


where UD is the detonation velocity of the explosive. The time at which the 
liner element reaches the axis is obtained from 


¢ (x) = T(x) +x tan a/V cos (a + 8) (46) 


and the point at which the liner element reaches the axis is given by 


z*(x) =x + (+7 - T)V sin (a + 8). (47) 


Values for the angle 5 can be computed from 


S = sin™*(V_cos € /20y)> (48) 

where tan(90° = a + «) = (L- H +x)/(x tan aq). (49) 
The depth of penetration at a time t will be given by 

E(t) = (b= t) 


where S is the standoff between the base of the charge and the top of the target. 
The jet velocity, Vile); appearing in the penetration equation is the value for 
the jet element at rhe stagnation point! between the jet and target. The jet 
velocity may also be regarded as a function of x, the original position of the 
element within the liner. The function V .(x) can be computed from the non-steady 
theory of jet formation’. The penetration velocity obtained by differentiating 
(50) is 


U = aP/at = V, + [+ . a(v,) | (av ,/at) , (51) 
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Ate fh. 7 ene, a) ee ae oe ae “ ee & 
atv) = (t' +V, at /av , - dz"/av,). (52) 


The differential equation obtained by eliminating U between (51) and (4) is 2 
av /at = -yw,/(1 +7) [t - av,)] . (53) 
J J an ee et 

The virtual origin approximation is obtained by lettin 
1 (53) is a first order linear differential equation of the form dv ,/at = 
d can be solved by well known numerical procedures. 

The numerical solution to the differential equation (53) requires a precise 
knowledge of the cone collapse parameters. The data presently available on the 


105mm unconfined test charge, while adequately describing the jet formation, do 


standoff equation. In lieu of using the generalization, the virtual origin 
approximation was shown to be satisfactory for standoffs down to half of a cone 
interest, and will 


be resolved when more precise data becomes available. 
EFFECTS OF COMPRESSIBILITY 


As pointed out in the introduction, the penetration of a target by a 


compressibilities of the target and jet are certainly not negligible at the 


pressures produced during the penetration, and it is important to estimate 


oa 


point with subsonic velocities. For adiabatic compressible flow, Bernoulli's 


equation requires that 


r Ps o° 0.4 
| Ps ap - (1/2)pcU = 0, (5h) 
O Q, 
rt 
rPs io 
ps 2 = 
= dp - (1/2) e-(V, - U) =0 (55) 
4 "9 v~ 
te) 
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aw 


: ° 
along the axial streamline, where p is the stagnation pressure; and o , the 


s ae 
respective densities at normal pressure. The effect of compressibiiity on the 
penetration velocity is readily determined by setting the left-hand side of (54) 


equal to the left-hand side of (55) and ere terms. Thus one obtains 


(1/2)o, (WV, a0) cote { [ (6° /;) (5/4) | dp. (56) 


If the material is incompressible the integral in (56) is zero and the resulting 
expression is equivalent to (1) with } set equal to one. It is apparent that the 
integral is zero if the jet and target are the same material. On the other hand, 
the largest error will result from treating as incompressible two materials that 


have greatly different compressibilities. 


It is possible to compute the error for any combination of materials for 

which suitable equation-of-state data are available. Usually, the high pressure 
dynamic measurements provide an empirical relation between the pressure and density 
along the shock Hugoniot. When the Griineisen ratio is also known, the pressure 

and density can be calculated for any adiabat that intersects the Hugoniot”. 

The penetration of a titanium target by a copper jet was chosen as an 

illustration because titanium is considerably more compressible than copper. 
Calculations were made using the equation of state data reported by McQueen and 
Marsh! Assuming a stagnation pressure of 500 kilobars, a penetration velocity 

of 0.412 em/u-sec is computed from (54), and a jet velocity of 0.726 cm/y-sec is 
obtained from (55). For a penetration ema of 0.412 em/-sec, the 
incompressible approximation predicts a jet velocity of 0.720 cm/p-sec, which 
is within 1% of the more exact calculations. The relative magnitudes of the 


quantities appearing in (56) are illustrated graphically in Figure 8. The 


For a general discussion concerning the behavior of solids at high dynamic 
pressures, the reader is referred to the section entitled "Compression of 
Solids by Strong Shock Waves" contributed by M. H. Rice, R. G. McQueen, and 
J. M. Walsh to Solid State Physics, Vol. 6, Edited by Fredrick Seitz and 


David Turnbull, Academic Press, Inc. Publishers (1958). 
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shaded area, is a small fraction of (1/2) ou, which is represented by the 
cross hatched area. It is noted, however, that the pressure computed from the 


incompressible approximation, 


aay 


(0/2) pur is 422 kilobars, ss than th 
actual pressure. Thus, as already demonstrated by direct experimental 
observations, the relationship between the jet and penetration velocities obtained 
is sufficiently accurate for all practical 


ate 4 7 CANNMAA GAY lM oO A. 4G 


problems involving the penetration of shaped-charge jets. 
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Line drawings illustrating the formation of a shaped- 
charge jet and the subsequent penetration of a target 
by the jet. Note the hydrodynamic nature of the 
penetration process. 
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Fig. 3 Penetration-time curve for the 105mm unconfined test charge fired at an 
8 1/4 in. standoff. Note that the continuous jet approximation is valid 
for the first 120u-sec. Beyond this point, the penetration time must 
be calculated using a method applicable to a jet of discrete particies. 
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Fig. 5 Curves representing the penetration produced by the 
br nee -up por ae ion of jets ee ae 105mm 
unconfined test charges fired at an 8 1/4 in. standoff. 
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Fig. 7 Line drawing illustrating the quantities employed in the 
hydrodynamic description of the collapse process for a 


lined conical cavity. 
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In earlier attempts to describe the penetration of shaped-charge jets, the 
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‘summation is taken over the penetration produced “by the discrete particles. It is . 
jshown that numerical values of the penetration variables for the 105mm unconfined 
icharge are compatible with the modified penetration equation. Scaling relations 


‘pre introduced, and the penetration equations are given in size independent 


variables. Although it was assumed that the incompressible eparccinat ton provides 
an adequate relationship between the jet and penetration velocities for a 
continuous jet, it is later shown that the finite compressibilities of the jet and 
jtarget materials do not seriously alter this relation. The incompressible 
lapproximation does underestimate the pressure by an appreciable amount. 
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